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NUMERICAL  MODELING  OF  FRINGING  FIELDS  AND  THEIR  USE 
FOR  COMPLEX  PERMimVITY  MEASUREMENTS  AT  HIGH  FREQUENCIES 


INTRODUCTION 

The  dielectric  measurements  of  biological  materials  are  increasingly  being 
performed  using  network  analysers  and  open-ended  coaxial  probes,  a  technique  first 
introduced  by  Burdette  et.  al.  in  1980  [1].  Since  then  there  has  been  an  increasing 
interest  in  the  use  of  open-ended  coaxial  probes  for  dielectric  measurements  because  of  its 
suitability  to  a  variety  of  applications  in  material  science,  chemical,  biochemical  and 
biomedical  fidds. 

Measurements  are  made  by  placing  the  probe  in  contact  with  a  sample  and 
measuring  its  admittance  or  reflection  coefficient  using  a  network  analyser  or  equivalent 
instrumentation.  Theoretically,  the  most  important  consideration  is  the  derivation  of  the 
relationship  between  the  measured  admittance  or  reflection  coefficient  and  the  dielectric 
properties  of  the  sample.  A  number  of  simplifying  assumptions  are  made  during  the 
course  of  the  derivation  which  restrict  the  validity  of  the  models  to  a  certain  range  of 
parameters.  We  have  examined  the  theoretical  basis  for  the  derivation  of  such  models  and 
their  suitability  for  dielectric  measurements  on  lossy  materials  in  the  frequency  range  300 
kHz  to  20  GHz. 

Two  situations  were  considered  depending  on  the  size  of  the  sample  in  contact  with 
the  prd)e.  First,  it  was  assumed  that  the  sample  is  homogeneous  and  occupies  the  entire 
half-space  beyond  the  probe.  The  second  case  considered  was  that  of  a  sample  of  finite 
thickness,  placed  between  the  probe  and  a  conducting  plane  parallel  to  it. 

The  penetration  of  the  electric  field  in  the  sample  was  determined  and  its 
dependence  on  the  field  frequency  and  the  dielectric  properties  of  the  sample  highlighted. 
The  study  of  the  field  penetration  is  important  because  it  helps  delineate  the  extent  of  the 
sampling  volume  for  a  particular  probe  and  frequency  range.  It  is  also  important  in  its 
own  right  and  may  be  used  to  assess  the  suitability  of  using  open-ended  coaxial  probes  for 
applications  other  than  dielectric  measurements. 


SEMI-INFINITE  SAMPLE 
Description  Of  The  Problem 

For  the  purpose  of  this  analysis  the  coaxial  probe  is  assumed  to  be  axisymmetric, 
fitted  with  a  ground  plane  that  is  large  by  comparison  to  its  diameter  and  that  in  addition 
to  the  principal  TEM  mode  only  TMon  evanescent  modes  that  preserve  cylindrical 
symmetry  are  generated  at  the  interface.  The  problem  is  best  described  in  the  cylindrical 
coordinate  system  (p,  (j>,  z)  as  illustrated  in  Figure  1  in  which  e  and  n  are  the  absolute 
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permittivity  and  permeability,  the  suffixes  0,  1  and  2  refer  to  air,  the  medium  inside  and 
outside  the  coaxial  line  respectively,  p  is  the  radial  distance  from  the  axis,  ((>  the  angular 
displacement  around  the  axis  and  z  the  displacement  along  the  axis.  In  the  first  instance 
we  assume  that  the  medium  outside  the  probe  is  uniform  and  occupies  the  entire  half¬ 
space  beyond  the  ground  plate. 

Expression  For  The  Hectric  And  Magnetic  Fields  Inside  The  Line 

The  fields  inside  the  line  are  a  superposition  of  the  forward  travelling  fundamental 
TEM  mode,  its  reflection  at  the  interface  and  the  TMot,  evanescent  modes  generated  at  the 
discontinuity.  The  field  components  can  be  derived  from  electromagnetic  theory  in  terms 
of  the  radial  potential  ^(p)  as 


hiP 


dp 


(1) 


Jl-1 


dp 


(2) 


n-l  '  In  / 


(3) 


in  which  F  is  the  reflection  coefficient  of  the  TEM  mode,  y„  is  the  propagation  constant 
of  mode  n 


and  k,^  =  ©Vi®!  is  the  wave  number,  ©  is  the  angular  frequency  and  q,  =  is 

the  intrinsic  impedance  of  the  line.  For  n  =  0  (TCM  mode)  Po  =  0,  for  n  >  0  p„  is  the 
nth  root  of 


Jo(P„b)Yo(p.a)-Jo(p„a)Yo(p„b)  =  0  (5) 

where  Jq  and  Yo  are  the  Bessel's  functions  of  the  first  and  second  kinds  of  order  zero,  a 
and  b  the  radii  of  the  coaxial  line. 

For  each  n,  the  function  is  obtained  from 

<I>n(P)  =  Jo(PnP)  Yo(p„a)-Jo(p„a)  Yo(p„p)  (6) 

Finally,  ao  and  a„  are  amplitude  terms  normalized  to  the  dimensions  of  the  line  and 
can  be  expressed  in  terms  of  the  magnitude  of  the  electric  field  at  the  aperture  E*p(p)  at 
z  =  0  as 
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In  this  derivation  [2-4],  the  open  end  of  the  coaxial  line  is  assumed  to  behave  as  a 
magnetic  source  radiating  into  the  z  >  0  region.  The  primed  coordinates  p*  and  <i>'  refer 
to  a  source  point  at  the  aperture  while  R  is  the  distance  from  the  source  to  the 
corresponding  field  point  which  at  z  =  0  is  R*  =  -1-  p’^  -  2pp'cos(<j)'-<|>). 


The  aperture  field  E*p(p)  at  z  =  0  can  be  obtained  from  (2)  as 


b;  ( p ) . -fililEL .  j; 


The  boundary  conditions  on  the  z  =  0  plane  require  that  the  tangential  electric  and 
magnetic  fields  to  be  continuous.  The  continuity  of  the  electric  field  is  ensured  by 
expressing  the  fields  on  either  side  in  terms  of  E*p.  From  (1)  and  (9),  the  continuity  of 
gives  after  substituting  for  Uq  and  F 


.  ^29 


in  which  Y*  is  the  terminating  admittance  of  a  coaxial  line  given  by  Y*  =  Yo(l-r)/(i+r) 
and  Yq  =  27r/[ii,ln(b/a)].  This  expression  was  first  derived  and  proposed  for  dielectric 
measurements  by  Mosig  et  al  (1981)  [S]. 

For  n  =  0,  TEM  propagation  is  assumed  throughout,  E*p(p)  is  inversely 
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proportional  to  p  and  expression  (11)  becomes 


y«=. 


i>l>  2s 


-jk^R 


[ln(jb/a)]2 


///  cos  ^  ^  d^'dp'dp 


«  «  0 


(12) 


By  expanding  the  exponential  term  in  (12)  into  an  infinite  series  and  dropping  the 
higher  order  terms  in  co  expression  (12)  becomes 


Y‘=jo>Cej  (13) 

in  which  C  is  a  frequency  independent  constant  whose  value  depends  on  the  dimensions  of 
the  probe.  By  virtue  of  its  derivation  (13)  is  strictly  valid  at  the  limit  of  low  frequencies. 
Expressions  (11)  and  (12),  derived  in  accordance  with  transmission  line  theory,  should  be 
expected  to  describe  more  adequately  the  aperture  admittance  of  the  probe. 

Expressions  (11),  (12)  and  (13)  are  collectively  described  as  models  for  the 
admittance  of  the  probe  in  terms  of  the  field  frequency  and  the  permittivity  of  the  material 
in  contact  with  it.  However,  only  the  simplest  formulation  (13)  has  an  analytical  solution, 
both  (11)  and  (12)  require  numerical  solutions. 


NUMERICAL  SOLUTIONS 

Numerical  solutions  are  approximate  techniques  that  can  reduce  the  rigor  of  the 
original  derivation  or  impose  numerical  limits  to  its  range  of  applicablility.  It  is  therefore 
important  to  report  on  the  numerical  implementation  for  models  (11)  and  (12). 

Method  Of  Moment  Solution 

The  exact  solution  of  (11)  requires  consideration  of  an  infinite  series  of  TMo„ 
modes  and  continuity  of  Uie  adds  over  the  entire  a  <  p  <  b  range.  In  practice,  only  a 
limited  number  of  modes  can  be  considered  and  matching  can  be  explicitly  satisfied  at  a 
limited  number  K  of  points  between  a  and  b  [5,6].  Because  of  the  axial  symmetry  this  is 
equivalent  to  matching  on  K  circles  within  the  aperture.  The  interval  [a,b]  is  thus  divided 
into  K  subintervals  with  boundaries  at  Pk  with  the  limits  Po  =  a  and  Pk  =  b.  The 
potential  in  each  ring  is  assumed  constant  with  a  value  such  that  pE*p(p)  =  V^.  Given 
that 

b 

T=-l*fE;{p)dpjv  ■  (14) 


where  V  is  the  incident  voltage  amplitude.  Expression  (1 1)  can  be  evaluated  at  the  mid¬ 
point  p,  =  (Pk+Pk.i)/2  of  every  K  interval  from 
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2V  ln(pt/pt.,) 

Hiln(l>/a)  i},ln(h/a) 


J^*a 

2% 


b2% 


j  j  cos  (♦''-♦) 


(15) 


Equation  (IS)  yields  a  set  of  K  linear  equations  in  V,^ 


2V  y  ^ 

qiln{i)/a)  ■**  * 


(16) 


that  can  be  expressed  in  the  matrix  form 


ra  =  IYJ(VJ  (17) 

and  solved  in  the  usual  way  for  V^.  The  reflection  coefficient  can  then  be  calculated  from 

r=-l*ig  Vn(p^/p,.j)  (18) 

In  principle  there  is  no  restriction  on  the  number  of  modes  that  are  incorporated  in 
the  solution,  in  practice,  including  more  than  10  modes  makes  almost  no  difference  to  the 
solution.  The  solution  is  more  sensitive  to  the  number  of  rings  K.  The  solution 
corresponding  to  K  =  oo  is  obtained  by  quadratic  extrapolation  [7]  from  solutions  at  3 
finite  values  of  K. 


Variational  Solutions 

There  are  three  possible  solutions  for  (12).  First,  the  approach  adopted  by  Misra 
(1987)  [8]  that  consists  in  expanding  the  exponential  term  into  a  series  and  performing  the 
triple  integration  on  the  individual  terms.  The  expression  for  the  admittance  of  the  probe 
takes  the  form 


n-l 


(19) 


where  the  terms  are  complex  terms  that  depend  on  the  dimensions  of  the  probe. 

On  evaluation  the  C2  term  goes  to  zero.  Misra  considered  only  the  first  three 
terms,  a  solution  that  proved  to  be  only  marginally  superior  to  the  simple  model  (13), 
However  this  technique  is  computationally  simple  and  more  terms  can  easily  be  added 
without  significant  time  penalty.  By  systematically  adding  more  terms  we  found  that  for 
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frequencies  in  the  GHz  region  and  lossy  materials  the  solution  was  sensitive  to  the  number 
of  terms  N  up  to  N  «  10. 

An  alternative  solution  of  ( '.2)  would  be  to  express  the  integral  with  respect  to  ()>' 
in  terms  of  Bessel  functions,  u*  uig  an  integral  representation  and  the  Neumann  formula 
[4],  which  gives  after  chan^mg  the  integration  parameters 


y-= 


j  YoJCa  r  [  Jq  ($jb)  -Jp  (^a)  ]  ^ 

ln{b/a)  J 


(20) 


where 

=  (0,  (21) 

I  n/2,  <  kj* 

The  integral  is  considered  separately  over  the  ranges  0  <  ^  <  kj  and 
k]  <  ^  <  00.  The  integrand  changes  from  imaginary  to  real  at  ^  =  kj,  and  the 
integration  over  these  two  ranges  gives  the  aperture  conductance  G  and  susceptance  B 
respectively.  For  0  <  ^  <  kj,  the  variable  of  integration  is  changed  to  ^  =  k2sin0  while 
for  k2  <  ^  <  00,  the  Dixon-Ferrar  integral  representation  for  the  product  of  two  Bessel 
functions  is  allied  to  transform  the  infmite  integral  into  a  finite  integral  of  integral  sines 
(Si).  It  follows  that 

Y*  «  G  +  jB  (22) 

where 

Yg  k; [ J-Q (k^asind)  - Jq (Jc^bsind)  ] ^ ^  ,23^ 

ln(jb/a)  /Cj  i  sin0 


(24) 


where  u*  =  a^  +  b^  -  2abcos4>  and  Si  is  the  sine  integral.  These  expressions  are  referred 
to  as  Macuvitz's  formulation.  Both  (23)  and  (24)  can  be  written  as  series  expansions  [9] 
that  are  convenient  for  numerical  calculation. 


The  third  and  most  straight  forward  solution  for  (12)  is  a  direct  evaluation  after 
performing  an  integration  transforma.  )n  to  convert  the  triple  integral  to  a  single  integral 
that  can  be  evaluated  using  Gaussian  integration  technique  [10].  This  is  the  technique  we 
advocate.  In  principle,  it  is  equivalent  to  the  other  techniques  when  the  number  of  terms 
of  their  respective  series  expansion  N  tmds  to  infinity.  To  illustrate  this  statement  the 
three  variational  techniques  were  used  to  calculate  the  admittance  of  a  Teflon  filled  probe 
of  radial  dimensions  a  =  0.456  mm,  b  =  1.490  mm,  in  contact  with  water,  in  the 
frequency  range  130  MHz  to  20  GHz.  The  results  showed  that,  for  N  =  10,  Marcuvitz’s 
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formulation  agreed  well  with  the  direct  solution;  the  exponential  expansion  solution  (19) 
agreed  with  them  for  frequoicies  up  to  7  GHz  beyond  which  it  bec^e  unstable. 


FINITE  THICKNESS  SAMPLE 

Expressions  (11)  and  (12)  can  be  adapted  to  the  case  of  a  dielectric  sample  of 
uniform  but  finite  thickness  backed  by  a  conducting  plane  parallel  to  the  ground  plane  of 
the  probe.  Both  planes  are  assumed  to  be  infinite  in  extent  which  is  equivalent  to 
neglecting  the  contribution  to  the  reflection  coefficient  of  the  field  behaviour  at  the  bound. 
One  way  of  performing  this  analysis  is  to  use  a  technique  known  as  the  method  of  images. 
The  sq)erture  of  the  probe  acts  as  a  source  of  magnetic  currents  that  gives  rise  to  fields  in 
the  region  of  the  didectric  sample  and  at  the  z  =  h  plane.  To  maintain  the  boundary 
conditions  at  z  =  0  and  z  -  h  require  an  infinite  sequence  of  images  of  the  source  at 
z  =  0  to  be  introduced  at  z  =  ±2h,  ±4h....  In  practice,  only  a  limited  number  of 
images  can  be  taken  into  consideration.  For  2q  images,  expressions  (11)  and  (12)  become 


e»,(p) 
Yx,  dp 


b2% 

a„+2ju€jj’  j E\{p')  cos  (♦'-♦)  p'^^'dp' 

s  0 


(25) 


bb2% 


with 


47li? 


[21) 


and  R,  is  the  distance  from  the  source  to  a  corresponding  field  point  which  at  z  =  2qh  is 

R,^  =  +  4q^h^. 


Expressions  (25)  and  (26)  can  be  solved  for  the  admittance  of  the  probe  using  the 
numerical  techniques  detailed  in  the  previous  section.  In  the  calculations  reported  here 
q  =  30,  the  solution  was  not  sensitive  to  further  increase  in  the  number  of  images. 


FIELD  PENETRATION  AND  ENERGY  ABSORPTION  IN  THE  SAMPLE 

It  is  evident  from  the  derivations  of  the  previous  sections  that,  except  at  the  limit 
of  low  frequencies,  the  field  in  the  z  >  0  region  is  due  to  radiation  from  the  principal 
TEM  and  higher  order  modes.  The  penetration  of  the  field  will  depend  on  the  frequency 
and  the  dielectric  properties  of  the  medium. 

The  field  is  derived  from  the  assumption  that  the  open  end  is  equivalent  to 
a  source  of  magnetic  current.  The  field  distribution  in  the  aperture  plane  takes 
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into  account  TEM  and  higher  order  modes.  For  a  practical  solution,  the  region 
between  inner  and  outer  conductors  is  subdivided  into  N  concentric  rings  of 

boimdaries  p^.i  and  Pn,/  there  is  a  constant  potential  within  each  ring.  The 
field  components  can  he  expressed  as 


1  A 

J  \za+jk,R)cos(i,'-^d(f,'dp' 


(28) 

(29) 


The  power  dissipated  in  the  sample  can  then  be  calculated  from 

P  =  a\E\^  (30) 


where  E  is  the  resultant  electric  field  and  a  the  total  conductivity  of  the  medium. 


RESULTS 

Semi-Infinite  Sample 

A  study  of  all  three  admittance  models  given  by  expressions  (11),  (12)  and 
(13)  revealed  the  systematic  difference  between  the  predictions  of  (11)  and  (12) 
due  to  the  contribution  of  higher  order  modes.  It  also  showed,  as  expected,  that 
the  simple  expression  for  the  admittance  (13)  dgrees  with  the  other  two  models 
only  at  low  frequencies  and  for  low  loss  low  permittivity  materials.  Model  (13) 
was  not  considered  further. 

The  predictions  of  models  (11)  and  (12)  were  compared  to  corresponding 
experimental  values  of  the  admittance  of  the  probe  in  contact  with  a  number  of 
standard  liquids.  The  measurements  were  made  using  a  network  analyser 
calibrated  using  standard  components  prior  to  connecting  the  probe.  Under  these 
conditions  the  experimental  values  are  expected  to  have  scatter  due  to  unwanted 
reflections  from  uncalibrated  components.  Typical  sets  of  results  are  given  in 
Figures  2  and  3.  In  Figure  2  the  real  and  imaginary  parts  of  the  admittance  of  the 
probe  in  methanol  are  plotted  as  a  function  of  frequency.  In  Figure  3  the  results 
for  water  and  formamide  are  given  in  the  complex  admittance  plane.  These 
results  show  that  the  experimental  values  do  not  agree  closely  with  either  model 
which  leads  to  the  conclusion  that  neither  model  is  suitable  for  absolute 
measurements.  However,  relative  measurements  can  be  simply  implemented  by 
performing  a  calibration  at  the  plane  of  the  probe  using  the  admittance  model  of 
choice  to  calculate  the  reflection  coefficient  of  the  standards  used.  This  method 
has  the  added  advantage  of  eliminating  all  unwanted  reflections.  Measurements 
using  this  technique  and  either  of  the  two  models  produces  accurate  dielectric 
measurements  in  the  frequency  range  up  to  20  GHz  as  evident  from  Figure  4 
which  shows  the  dielectric  spectrum  of  two  standard  liquids  measured  using  (11) 
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and  (12).  An  analysis  of  these  data  showed  that  the  two  set  measurements  are 
not  significantly  different  from  each  other  and  from  the  previously  published 
and  accepted  values. 


Models  (11)  and  (12)  are  therefore  equally  suitable  for  performing  dielectric 
measurements  at  frequencies  up  to  20  GHz.  The  choice  between  them  may 
therefore  be  made  on  the  grounds  of  convenience  since  neither  has  been  shown 
to  be  superior  to  the  other.  Model  (11)  requires  a  rather  complex  numerical 
solution  for  its  implementation  and,  for  this  reason,  is  unsuitable  for  real  time 
measurements.  By  contrast,  model  (12)  can  be  solved  directly  using  integral 
transformation  and  Gaussian  integration  techniques.  To  show  the  potential  of 
the  technique  measurements  were  made  on  electrolyte  solutions  ranging  in 
concentration  from  very  dilute  to  the  saturated  solution.  The  results  for  KCl  at 
20  *C  are  given  in  Table  1  as  parameters  of  the  Cole-Cole  relationship 


£  =e  + 


1-a 


\  +  {jQ}Z)  G)£, 


(31) 


in  which  o  is  the  ionic  conductivity  and  all  the  parameters  have  their  usual 
significance.  The  results  for  NaCl  solutions  are  given  in  Table  2.  The  study  of 
the  variation  of  the  parameters  with  concentration,  which  is  not  simply 
monotonic,  carries  information  on  the  structure  of  the  water  in  the  solutions 
and  on  the  presence  and  nature  of  water  bound  to  the  electrolyte.  Examples  of 
the  concentration  dependence  of  the  dielectric  parameters  are  as  follows: 

For  NaCl  at  25 ’C 


e,  =  78. 3  - 1. 86x -t- 2. 84  X 1 0'^ 
t  =  8.1  +  7.26x10-^;c-1.04x10'';c'-h6.14x10'^x' 
(T=1.61jc-2.64x10-^a:^ 
a  =  8.68  X  IQ-^x  +  2.39  x 

For  KCl  at  20 ’C 

£,  =  80.1  -  1.60:c-h  4.18  X 10'^  x' 

T  =  9.3  -  5.69  X 10'^  jc  -  6.84  X  10'';c^  +  6.03  x  lO^^ 
ct=1.72x-1.46x10'^x^ 
a  =  5. 39  X 10'^  X  +  2.5 1 X 10“^ 

where  x  is  the  concentration  expressed  as  percentage  by  weight,  the  relaxation 
time  is  in  pico-seconds  and  the  conductivity  in  siemens  per  metre.  Previous  data 
on  the  dielectric  properties  of  ionic  solution  was  available  in  a  1971  publication 
by  Stogryn  [10].  Stogryn's  models  were  based  on  limited  data  and  a  Debye  model 
was  assumed  without  justification.  The  analysis  of  our  data  is  still  underway. 
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Finite  Thickness  Sample 


Measurements  performed  using  the  finite  thickness  technique  showed 
that  for  sample  thicknesses  h  ^  (b  -  a)  results  and  conclusions  similar  to  those 
obtained  for  the  semi-infinite  sample  can  then  be  reached.  Both  (25)  and  (26) 
have  been  used  successfully  to  perform  dielectric  measurements  and,  provided 
that  h  ^  (b  -  a),  the  accuracy  of  the  measurements  is  not  compromised  by  the 
choice  of  model.  Table  3  gives  the  dielectric  parameters  of  4  standard  liquids 
measured  using  this  technique.  The  technique  reverts  to  the  case  of  a  semi 
infinite  sample  when  h  tends  to  infinity.  The  dielectric  parameters  in  Table  2 
correspond  to  those  of  the  general  relationship 


£  —£ 


[l  +  (ya)T)'-“f  (»£o 


(34) 


which  correspond  to  the  Debye  equation  for  a  =  0  and  P  =  1.  It  corresponds 
to  the  Cole-Cole  equation  for  P=1  and  0<a<l  and  to  the  Cole-Davidson 
relationship  for  a  =  0  and  0<P<1. 

Field  Penetration  And  Power  Dissipation 

Calculations  were  made  using  (28),  (29)  and  (30)  for  a  Teflon  filled  probe  of 
radial  dimensions  a  =  0.456  mm,  b  =  1.490  mm  for  methanol,  water  and  0.9% 
NaCl  solution  (physiological  saline)  at  three  frequencies  200  MHz,  2.45  GHz  and 
20  GHz.  The  results  are  presented  as  follows: 

Figure  5  a,  b  and  c  show  the  normalised  power  dissipation  as  a  function  of 
p,  at  z  =  0.1mm,  in  the  three  samples,  at  the  three  frequencies  respectively. 

Figure  6  a,  b  and  c  show  the  normalised  value  of  |£|^as  a  function  of  p,  at  z 
=  0.1mm,  at  the  three  frequencies,  in  the  three  samples  respectively. 


Figure  7  a,  b,  and  c  show  the  normalised  value  of  |£f ,  as  a  function  of  p, 
for  z  values  ranging  from  0.1  to  3mm,  at  the  three  frequencies  respectively. 

These  results  show  that  in  all  cases,  in  the  radial  direction,  the  maximum 
power  absorption  occurs  at  the  edge  of  the  centre  conductor.  Beyond  the  centre 
conductor,  the  power  dissipated  decreases  exponentially  to  a  constant  value  as  it 
approaches  the  radius  of  the  outer  conductor  b,  beyond  which  it  suffers  another 
exponential  decay  to  vanishingly  small  values. 

At  a  given  distance  (z=0.1mm)  from  the  probe  the  power  deposition 
depends  on  the  dielectric  properties  of  the  sample.  For  example,  at  200  MHz  the 
highest  power  absorption  occurs  in  NaCl  due  to  its  high  conductivity.  At  20  GHz 
the  highest  absorption  occurs  in  methanol  due  to  high  field  penetration. 
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From  Figures  7  it  can  be  seen  that  the  bulk  of  power  is  absorbed  in  the 
terminating  lossy  media  within  the  hemisphere  subtended  by  the  probe.  Outside 
this  volume  the  power  absorbed  in  the  samples  is  less  than  -25  dB  at  200  MHz 
and  2.45  GHz  but  less  than  -20  dB  at  20  GHz.  This  indicates  that  the  ground  plate 
is  more  important  at  the  highest  frequency. 


CONCLUSIONS 

The  main  purpose  of  this  project  was  to  develop  an  accurate  technique  to 
measure  the  complex  permittivity  of  materials  based  on  a  rigorous  theoretical 
formulation  of  the  admittance  of  an  open-ended  coaxial  probe.  This  has  been 
achieved.  The  technique  was  tested  by  measurements  on  standard  solution. 

The  new  technique  was  used  to  perform  a  comprehensive  study  of  ionic 
solutions.  The  extensive  data  obtained  led  to  the  development  of  models  to 
predict  the  dielectric  parameters  of  ionic  solutions  at  a  given  temperature  as  a 
function  of  concentration. 

The  theoretical  study  and  numerical  techniques  were  extended  to  treat  the 
case  of  a  thin  sample  backed  by  a  metallic  plane.  The  technique  was  tested  with 
standard  solutions. 

Further  analysis  led  to  the  development  of  expressions  for  the  field  in  the 
sample  adjacent  to  the  probe.  This  enabled  the  pattern  of  power  deposition  to  be 
determined  and  improved  our  understanding  of  the  sampled  volume.  More 
work  is  needed  to  determine  the  power  deposition  in  biological  material. 
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Figure  3.  Complex  admittance  of  probe  in  (a)  water,  (b)  formamide. 
(The  numbers  are  frequencies  in  GHz) 
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TABLE  1.  Dielectric  parameters  of  KCl  solutions  obtained  by  analysis  of  the 

experimental  results  at  20 *C.  The  A  terms  correspond  to  the  95% 
confidence  interval. 


Cone.  % 

ES 

Aes 

x(ps) 

Ax  (ps) 

a 

Aa 

a  (Sm-i) 

Aa  (Sm-i) 

25.6 

66.83 

1.45 

0.306 

0.020 

34.528 

0.067 

20.5 

64.48 

0.27 

10.30 

0.11 

0.218 

0.006 

29.271 

0.015 

15.3 

65.34 

0.33 

8.89 

0.13 

0.142 

0.010 

22.967 

0.017 

1Z8 

66.21 

0.25 

8.80 

0.10 

0.096 

0.008 

19.576 

0.011 

10.2 

68.09 

0.20 

8.75 

0.08 

0.071 

0.007 

16.021 

0.009 

7.67 

70.53 

0.16 

8.77 

0.06 

0.054 

0.005 

12.250 

0.007 

5.12 

73.30 

0.23 

8.83 

0.09 

0.045 

0.007 

8.447 

0.008 

2.56 

76.44 

0.11 

8.99 

3.04 

0.031 

0.003 

4.456 

0.004 

2.05 

77.08 

0.09 

9.09 

0.03 

0.017 

0.003 

3.552 

0.003 

1.53 

77.81 

0.08 

9.16 

0.03 

0.014 

0.002 

2.712 

0.003 

1.28 

78.13 

0.04 

9.06 

0.02 

0.023 

0.001 

2.334 

0.001 

1.02 

78.67 

0.13 

9.18 

0.05 

0.016 

0.004 

1.849 

0.003 

0.767 

79.15 

0.12 

9.19 

0.05 

0.021 

0.004 

1.418 

0.002 

0.512 

79.43 

0.08 

9.20 

0.04 

0.012 

0.003 

0.959 

0.001 

0.256 

79.81 

0.08 

9.21 

0.04 

0.012 

0.003 

0.500 

0.002 

0.205 

80.00 

0.07 

9.29 

0.03 

0.009 

0.002 

0.404 

0.001 

0.153 

80.01 

0.06 

9.26 

0.03 

0.011 

0.002 

0.307 

0.001 

0.128 

79.94 

0.06 

9.19 

0.02 

0.013 

0.002 

0.258 

0.001 

0.102 

80.09 

0.07 

9.26 

0.03 

0.012 

0.002 

0.208 

0.001 

0.0767 

79.91 

0.05 

9.15 

0.02 

0.011 

0.001 

0.159 

0.001 

0.0512 

80.02 

0.06 

9.21 

0.02 

0.008 

0.002 

0.107 

0.001 

0.0256 

80.01 

0.05 

9.17 

0.02 

0.009 

0.001 

0.054 

0.001 

TABLE  2.  Dielectric  parameters  of  NaCl  solutions  obtained  by  analysis  of  the 

experimental  results  at  25 *C.  The  A  terms  correspond  to  the  95% 
confidence  interval. 


Cone.  % 

es 

Aes 

x(ps) 

Ax  (ps) 

a 

Aa 

a  (Sm'b 

Aa  (Sm-i) 

24.65 

1.37 

11.20 

0.93 

0.354 

0.028 

23.521 

0.053 

20.6 

52.51 

0.93 

9.10 

0.44 

0.287 

0.023 

21.976 

0.046 

10.89 

61.21 

0.27 

7.79 

0.10 

0.112 

0.009 

14.277 

0.016 

5.64 

68.52 

0.15 

7.92 

0.05 

0.048 

0.005 

8.166 

0.009 

3.77 

72.08 

0.21 

8.00 

0.09 

0.052 

0.007 

5.701 

0.007 

1.91 

74.90 

0.13 

8.05 

0.05 

0.021 

0.004 

3.085 

0.004 

0.73 

77.07 

0.12 

8.04 

0.05 

0.017 

0.004 

1.266 

0.002 
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TABLE  3.  Dielectric  parameters  of  standard  liquids  obtained  by  analysis  of  the 

experimental  results  on  finite  thickness  samples  at  20  *C.  The  A  terms 
correspond  to  the  95%  confidence  interval. 


h(mm) 

Es 

AEs 

e. 

Ae«, 

t(ps) 

AT(ps) 

WSm'b 

AO(Sin-b 

AP 

Methanol 

1.0 

33.8 

0.2 

5.0 

0.3 

50.2 

1.1 

- 

- 

- 

- 

1.5 

35.3 

0.2 

4.6 

0.2 

54.0 

0.7 

- 

- 

- 

- 

2.0 

34.6 

0.1 

4.4 

0.1 

53.3 

0.4 

- 

- 

- 

- 

3.0 

34.7 

0.1 

4.6 

0.1 

54.0 

0.3 

- 

- 

- 

- 

oo 

33.7 

0.03 

4.85 

0.04 

53.8 

0.2 

- 

- 

- 

- 

Ethanol 

1.0 

25.3 

0.3 

4.6 

0.2 

173 

5 

1.5 

27.2 

0.2 

4.4 

0.2 

164 

4 

- 

- 

- 

- 

2.0 

25.7 

0.2 

4.2 

0.1 

174 

3 

- 

- 

- 

- 

3.0 

26.2 

0.2 

4.2 

0.1 

176 

3 

- 

- 

- 

- 

eo 

25.1 

0.1 

4.22 

0.08 

179 

2 

- 

- 

- 

- 

Ethanediol 

1.0 

42.4 

0.2 

3.1 

0.3 

170 

5 

- 

- 

0.73 

0.02 

1.5 

43.3 

0.2 

3.0 

0.3 

162 

5 

- 

- 

0.75 

0.02 

2.0 

42.3 

0.1 

3.4 

0.2 

156 

3 

- 

- 

0.77 

0.02 

3.0 

42.6 

0.1 

3.6 

0.2 

158 

3 

- 

- 

0.82 

0.01 

OO 

42.0 

0.01 

3.66 

0.02 

162 

0.3 

- 

- 

0.806 

0.001 

Formamide 

1.0 

110.3 

0.4 

5.2 

0.7 

36.8 

0.5 

0.037 

0.007 

- 

- 

1.5 

112.1 

0.5 

5.3 

0.8 

38.1 

0.6 

0.038 

0.008 

- 

- 

2.0 

110.8 

0.3 

5.5 

0.4 

36.9 

0.3 

0.035 

0.004 

- 

- 

3.0 

111.4 

0.2 

5.1 

0.3 

37.7 

0.2 

0.031 

0.003 

- 

- 

OO 

111 

0.06 

5.7 

0.1 

40.0 

0.08 

0.028 

0.001 

- 

- 
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